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ABSTRACT 

Study of the neutrinoless double beta decay, Qv(3fi, includes a variety of problems of 
nuclear structure theory. They are reviewed here. The problems range from the mech- 
anism of the decay, i.e. exchange of the light Majorana neutrino neutrino vesus the 
exchange of some heavy, so far unobserved particle. Next, the proper expressions for 
the corresponding operator are described that should include the effects of the nucleon 
(*~^ 1 size and of the recoil order terms in the hadronic current. The issue of proper treament 

fSJ . of the short range correlations, in particular for the case of the heavy particle exchange, 

£jf^ is discussed also. The variety of methods employed these days in the theoretical evalua- 

^3 ■ tion of the nuclear matrix elements M 0v is briefly described and the difficulties causing 

' the spread and hence uncertainty in the values of M 0v are discussed. Finally, the issue 

of the axial current quenching, and of the resonance enhancement in the case of double 
electron capture are described. 

1. Introduction 

Study of the neutrinoless double beta decay is the most sensitive current test of the total lepton number 
conservation. If observed, it would serve as a proof that not only the lepton number Ltot is not conserved, 
but also that neutrinos are massive Majorana fermions. Hence the study of this process is one of the 
important parts of the search for "Physics Beyond the Standard Model" . Different aspects of the complex 
of problems associated with the (3(3 decay are described in this focus issue. Since the process involves 
1 nuclei, i.e. complicated many body systems, the analysis of the decay necessary involves various nuclear 
Q\, structure problems. In this article the general questions of nuclear structure, relevant for the understanding 
. of the (3(3 decay rate, are discussed. In other works in this issue the particular approximate techniques 
OO \ for the treatment of the corresponding nuclear matrix elements are described. Here I try to discuss the 
more general framework. Various aspects, both theoretical and experimen tal, of the (3(3 decay have been 
. reviewed many times. I quote here just some of the review articleJU2I2I3I3 ; earlier references can be found 
| there. 

In double beta decay two neutrons, bound in the ground state of an even-even initial (or parent) 
nucleus are transformed into two bound protons, typically also in the the ground state of the final (or 
| granddaughter) even-even nucleus, with the simultaneous emission of two electrons only for the 0^(3(3 
mode, or two electrons plus two v e for the 1v(3(3 mode. Transitions leading to the excited bound states 
of the final nucleus are sometimes kinematically allowed as well; however, we will concentrate here on the 
most often studied case of the ground state to ground state decays 

(Z, A) g . s . + A) g . s . + 2e" + (2P e ) . (1) 

Note that transitions of two bound protons into two bound neutrons (nuclear charge is decreased by 
two units in that case) with the emission of either two positrons, or a positron accompanied by the electron 
capture, or by two simultaneous electron captures, are also kinematically allowed for several even-even 
nuclei, again in both two-neutrino and neutrinoless modes. The lepton phase space for these transitions is 
suppressed compared to the decays of the type ([TJ, while the nuclear structure issues are analogous. Hence, 
again we will concentrate on the decays described in Eq. (pQ). However, the QvECEC mode in few cases 
can have a resonance character (Q — >■ 0). The possible enhancement in those special cases is described in 
the Appendix 

The (3(3 decay, in either mode, can proceed only if the initial nucleus is stable against the standard (3 
decay (both (3~ and (3 + or EC). That happens exclusively in even-even nuclei, where moreover the ground 
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state is always I w = + . The (3/3 decay rate is a steep function of the energy carried by the outgoing leptons 
(i.e. of the decay Q-value). Hence, transitions with larger Q-value are easier to observe. For this reason 
in Table El I list all candidate nuclei with Q values larger than 2 MeV that are particularly well suited for 
the study of the (3(3 decay. 

There has been a significant progress recently in the accuracy of the atomic mass determination using 
various trap arrangements. In many cases the Q-values are determined with accuracy better than 1 keV, 
making the search for the all important Qv(3(3 decay mode easier; in Table El these more recent Q-value 
determinations are shown, together with the corresponding references. 

In both modes of the (3(3 decay the rate can be expressed as a product of independent factors that 
depend on the atomic physics (the so called phase-space factors G 0u and G 2v ) that include also the Q-value 
dependence as well as the fundamental physics constants, nuclear structure (the nuclear matrix elements 
M 0l> and M 2v \ and for the 0v(3(3 mode the possible particle physics parameters (the effective neutrino 
mass {rapp) in the simplest case). Thus 

1 = G ^ ]M ^\ 2 \{m^)\ 2 ; -4r = G 2u \M 2u \ 2 . (2) 

1 l/2 1 l/2 



Table 1: Candidate nuclei for (3(3 decay with Q > 2 MeV 
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The values of G 0u and G 2v are also listed in Table [TJ The entries there are taken from Ref. ' , and 
are not corrected for the small changes in Q and gA since that time. Also, since by convention the nuclear 
matrix elements M 0u are dimensionless, the nuclear radius appears in them as a multiplicative factor. To 
compensate for it, the phase-space factor G 0u is proportional to R~ 2 , where R = x A 1 / 3 is the nuclear 
radius. In Table [1] the value tq = 1.2 fm was used. (Note that, obviously, the values of the phase-space 
factors depend on the convention used for ro and gA- One has to keep that issue in mind when using the 
Eq. to relate the half-lifes and nuclear matrix elements (see e.g. H^ 7 - 1 ).) 

Double beta transitions are possible and potentially observable because nuclei with even Z and N are 
more bound than the odd-odd nuclei with the same A = N + Z. A typical example is shown in Fig. [TJ 
With one exception, all nuclei in Table E] have an analogous mass pattern. The one exception is 48 Ca 
where the intermediate nucleus 48 Sc can be in principle reached by the (3~ decay of 48 Ca with Q= 278 
keV. However, the ground state of 48 Sc is 6 + and the first excited state at 131 keV is 5 + . (3 decays with a 
large nuclear spin change are heavily suppressed; in this particular case the (3~ decay of 48 Ca has not been 
observed as yet, while the 2v(3(3 decay has been observed. 

The two-neutrino mode (2v(3(3) is just an ordinary beta decay of two bound neutrons occurring simul- 
taneously since the sequential decays are forbidden by the energy conservation law. For this mode, clearly, 
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Figure 1: Atomic masses of the isotopes with A — 136. Nuclei Xe, Ba and Ce are stable against the ordinary /3 
decay; hence they exist in nature. However, energy conservation alone allows the transition 136 Xe — > 136 Ba + 2e~ (+ possibly 
other neutral light particles) and the analogous decay of 136 Ce with the positron emission. 



the lepton number is conserved and this decay is allowed in the standard model of electroweak interaction. 
It has been repeatedly observed in a number of cases and proceeds with a typical half-life of ~ 10 19-20 years 
for the nuclei with Q-values above 2 MeV. In contrast, the neutrinoless mode obviously violates the 

law of lepton number conservation and is forbidden in the standard model. Hence, its observation would, 
as already stated, be a signal of the " physics beyond the standard model". 

The two modes of the f3f3 decay have some common and some distinct features. The common features 

are: 

• The leptons carry essentially all available energy. The nuclear recoil is negligible, Q/Am p >C 1. 

• The transition involves the + ground state of the initial nucleus and (in almost all cases) the + 
ground state of the final nucleus. In few cases the transition to an excited + or 2 + state in the final 
nucleus is energetically possible, but suppressed by the smaller phase space available. (But the 2^/3/3 
decay to the excited + state has been observed in few cases.) 

• Both processes are of second order of weak interactions, ~ G F , hence inherently slow. The phase 
space consideration alone (for the 2uj3j3 mode ~ Q 11 and for the Ov/3/3 mode ~ Q 5 ) give preference to 
the Oz//3/3 which is, however, forbidden by the lepton number conservation and therefore much slower 
(at least by a factor 10 5 for the nuclei listed in Table [T]) than the 2v(3(3 decay. 

The distinct features are: 

• In the 2u/3f3 mode the two neutrons undergoing the transition are uncorrelated (but decay simulta- 
neously) while in the Ou/3/3 the two neutrons are correlated. In the 2vf3f3 mode the corresponding 
momentum transfer q is restricted by the decay Q- value; hence qR 1. On the other hand, in the 
Qvf3f3 mode the momentum transfer is of the order of the nucleon Fermi momentum q ~ qFermi'i 
hence qR > 1 in that case. 

• In the 2vf3f3 mode the sum electron kinetic energy T\ + T2 spectrum is continuous and peaked below 
Q/2. This is due to the electron masses and the Coulomb attraction. As T\ + T2 — > Q the spectrum 
approaches zero approximately like (A.E/Q) 6 . 
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• On the other hand in the Ovf3f3 mode the sum of the electron kinetic energies is fixed, T\ + T 2 = Q, 
smeared only by the detector resolution. This allows one to separate the two modes experimentally 
by measuring the sum energy of the emitted electrons with a good energy resolution, even if the 
decay rate for the mode is much smaller than for the 2vj3f3 mode. 

Another hypothetical mode of double-beta decay is often considered in the literature, the decay accom- 
panied by Majoron emission. Majorons are supposed to be very light or massless particles x° that couple 
to neutrinos. A variety of approaches involving massless (or almost massless) scalar particles and their 
impact on Ov/3/3 decay have been considered. The discussion of this topic goes beyond the scope of the 
present review; a rather complete list of references can be found e.g. in Ref. In all of these hypotheses 
the sum electron spectra are continuous with the characteriostic shape (neglecticng Coulomb effects) 

~ (Q - E 1 - 2 ) n p lP2 E 1 E 2 , (3) 



dE 1 dE 2 

where the "index" n = 1 — 7 depends on the model in question. The decay rate is expressed as 

T°» M = \(g x )\ 2 \M x \ 2 G x (Q,Z) , (4) 

where {g x ) is the model dependent averaged coupling constant. For n = 1 that constant is experimentally 
constrained to be less than about 10~ 5 . 



2. Mechanism of the 0^/3/3 decay 

The relation between the Oz^/3/3-decay rate and the effective Majorana mass {mnpp) is to some extent 
problematic. The rather conservative assumption leading to Eq.([2]) is that the only possible way the Oz^/3/3 
decay can occur is through the exchange of a virtual light, but massive, Majorana neutrino between the 
two nucleons undergoing the transition, and that these neutrinos interact by the standard left-handed 
weak currents. But that is not the only theoretically possible mechanism. Lepton number violating (LNV) 
interactions involving so far unobserved much heavier (~ TeV) particles might lead to a comparable Oz^/3/3 
decay rate. 

In general Oz^/3/3 decay can be gener ated by (i) light massive Majorana neutrino exchange or (ii) heavy 
particle exchange (see, e.g. resulting from LNV dynamics at some scale A above the elec- 

troweak one. The relative size of heavy (Ah) versus light particle (A£) exchange contributions to the decay 
amplitude can be crudely estimated as follows : 

L (k 2 ) ' H A5 ' A L A 5 (mpp) ' { 1 

where (map) is the effective neutrino Majorana mass, (k 2 ) ~ (100 MeV) 2 is the typical light neutrino 
virtuality, and A is the heavy scale relevant to the LNV dynamics. Therefore, Ah/Al ~ 0(1) for (mpp) ~ 
0.1 — 0.5 eV and A ~ 1 TeV, and thus the LNV dynamics at the TeV scale leads to similar Oi//3/3-decay 
rate as the exchange of light Majorana neutrinos with the effective mass (mpp) ~ 0.1 — 0.5 eV. 

Obviously, the Qv(3(3 lifetime measurement by itself does not provide the means for determining the 
underlying mechanism. The spin-flip and non-flip exchange can be, at least in principle, distinguished by the 
measurement of the single-electron spectra or polarization. However, in most cases the mechanism of light 
Majorana neutrino exchange, and of heavy particle exchange cannot be separated by the observation of the 
emitted electrons. We will not discuss here the possible ways of determining which mechanism is responsible 
for the transition, once the Oz^/3/3 decay has been observed. However, obviously, the corresponding nuclear 
matrix elements M 0u depend on that. 

Lets comment now on the main differences between the mechanism involving light or heavy particle 
exchange. We will show later that the evaluation of the nuclear matrix element M 0u can be performed in 
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the closure approximation, i.e. without explicit treatment of the virtual states in the intermediate odd-odd 
nucleus. Thus 

M°» = (f\\0 K \\i) , (6) 

where the operator O creates two protons and annihilate two neutrons. In addition, the operator O k 
depends on the distance between these nucleons, and on their other quantum numbers. The main difference 
between the mechanism involving the light massive Majorana neutrino exchange and the heavy (~ TeV) 
particle exchange is the range of the operator K . The light neutrino exchange represents two point-like 
vertices separated by the distance r ~ 1/q. The decay rate is then proportional to the square of the effective 
Majorana neutrino mass as in Eq. ([2]). On the other hand the heavy particle exchange represents a single 
point-like vertex (six fermions, four hadrons and two leptons), i.e. dimension 9 operator. Proper treatment 
of the short range nucleon-nucleon repulsion is obviously crucial in that case. The relation between the 
neutrino mass and the decay rate is not simple in that case. 

Independently of its mechanism the existence of 0u(3(3 decay would mean that on the elementary particle 
level a six fermion lepton number violating amplitude transforming two d quarks into two u quarks plus 
two electrons is nonvanishing. As was first pointed out by Schechter and Valle^^ 1 , already thirty years ago, 
this fact alone would guarantee that neutrinos are massive Majorana fermions. This qualitative statement 
(or theorem), however, as we pointed out above, does not in general allow one to deduce the magnitude 
of the neutrino mass once the rate of the Ov(3(3 decay have been determined. It is important to stress, 
however, that quite generally an observation of any total lepton number violating process, not only of the 
decay, would necessarily imply that neutrinos are massive Majorana fermions. 



3. 2vfif} decay 

Study of the 2v(3(3 decay is an important nuclear physics problem by itself. Moreover, evaluation of 
the M 2v matrix elements is an important test for the nuclear theory models that aim at the determination 
of the analogous but different quantities for the more fundamental Of neutrinoless mode. So, we begin our 
discussion with that experimentally more accessible (3(3 decay mode. 

The rate of the 2v(3(3 decay was first estimated by Maria Goeppert-Meyer already in 1937 in her thesis 
work suggested by E. Wigner, basically correctly. Yet, first experimental observation in a laboratory 
experiment was achieved only in 1987, fifty years later . (Earlier observations of the (3(3 decay E21 ? > ? ) 
were based on the geochemical method that cannot separate the 2u(3(3 and 0v(3(3 decay modes. Later 
laboratory measurements have shown, as expected, that the geochemical method determines dominantly 
the 2v(3(3 decay mode.) 

Note that such delay is not really exceptional in neutrino physics. It took more than twenty years 
since the original suggestion of Pauli to show that neutrinos are real particles in the pioneering experiment 
by Raines and Cowan. And it took another almost fifty years since that time to show that neutrinos are 
massive fermions. Why it took so long in the case of the (3(3 decay? As pointed out above, the typical 
half-life of the 2v(3(3 decay is ~ io 19-20 years. Yet, its "signature" is very similar to natural radioactivity, 
present to some extent everywhere, and governed by the half-life of ~ 10 10 years, or much less for most of 
the man-made or cosmogenic radioactivities. So, background suppression is the main problem to overcome 
when one wants to study either of the (3(3 decay modes. 

During the last two decades the 2v(3(3 decay has been observed in "live" laboratory experiments in 

many nuclei, often by different groups and using different methods. That shows not only the ingenuity 

of the experimentalists who were able to overcome the background nemesis, but makes it possible at the 

same time to extract the corresponding 2v nuclear matrix element from the measured decay rate. The 

resulting nuclear matrix elements M 2u , which have the dimension energy -1 , are plotted in Fig|2l Note the 

pronounced shell dependence; the matrix element for 100 Mo is almost ten times larger than the ones for 
i30 Te Qr i36 Xe _ 
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Figure 2: Matrix elements M 2v in MeV 1 based on the experimental halflife measurements. 



The derivation of the 2i//3/3-decay rate formula is analogous to the treatment of ordinary beta decay. 
It begins with the Fermi golden rule for second order weak decay 
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(7) 



where the sum over m includes all relevant virtual states in the intermediate odd-odd nucleus and fi labels 
the different Dirac structures of the weak interaction Hamiltonian. 

When calculating the + — > + transitions, it is generally a very good approximation to replace the 
lepton energies in the denominator by the corresponding average value, i.e., E e +p u ~ Eq/2, where the 
Eq = Mi — Mf is the total decay energy including electron masses. The lepton momenta, for both electrons 
and neutrinos are all q < Q and thus qR <C 1, where R is the nuclear radius. Hence the so-called long 
wavelength approximation is valid and the rate formula separates into a product of the nuclear and lepton 
parts, where the lepton part G 2u (Eq, Z) is just the phase-space integral discussed and tabulated earlier. 

The nuclear structure information is contained in the nuclear matrix element; only the Gamow- Teller 
or part contributes in the long wavelength approximation 
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The individual terms in the eq. ([8]) have a well defined meaning, in particular for the most relevant ground 
state to ground state transitions. The terms (m\ffk T k\^t) represent the amplitudes of the (3~ strength 
of the initial nucleus and can be explored in the nucleon charge exchange reactions such as (p, n) and 
( 3 He,t). On the other hand the terms (0~j\aiT^\m) represent the /3 + strength in the final nucleus and can 
be explored in the opposite nucleon charge exchange reactions such as (n,p) and (d, 2 He). In this way one 
can (up to the sign) explore the contribution of several low lying states to the M 2v matrix element. (See 
the article by H. Ejiri and D. Frekers in this issue.) 

In Eq.([8]) the energy denominators signify that the low- lying intermediate I w = 1 + states contribute 
significantly more that the high- lying states. Nevertheless, it is useful to analyze a related expression 
obtained in the closure approximation: 
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where AE is the average energy denominator defined by the above equation. The numerator of Eq.Q is 
the definition of the closure 2v nuclear matrix element. By itself, the 2v closure matrix element does not 
have a direct physics interpretation. 

For the two-body operator (i.e. in closure) in the 2v matrix element one can derive an approximate 
sum rule, analogous to the famous Ikeda sum rule of the GT operator . Remember that for the beta 
strengths = SffjTj the Ikeda sum rule is obeyed independently of the nuclear structure as long as 
nuclei are assumed to be made of nucleons only 

s p- _ 5 /3+ = 3^ _ Z ) . (10) 

In nuclei with neutron excess, N > Z the /3 + strength is much smaller than the (3~ strength, which is 
concentrated in the giant GT resonance. 

In analogy, the strength corresponding to the double beta operator ^i 1 k^i T i r ^k' T k wu ^ ^ e concentrated 
mostly in the "double GT resonance" and its strength will be approximately equal to the sum rule 6(iV — 
Z)(N — Z + 1)1251 . Assuming further that the average energy denominator AE in Eq. ([9|) is 0(1) MeV, we 
see that 2^/3/3 closure matrix element, connecting the ground states of the initial and final nuclei (see Fig. 
[2]) is heavily suppressed, representing only a very small fraction of the corresponding sum rule. From that 
it follows that an accurate evaluation of the matrix elements M 2u and M 2 ? is going to be quite difficult, 
since these quantities are so small in their natural units. 

To shed more light on the problem, consider the dependence of the 2^/3/3 matrix elements on the 
distance between the two neutrons that are transformed into two protons. These two neutrons are not 
correlated, but nevertheless they are both bound in the corresponding nuclei, and decay together. To 
characterize the radial dependence lets introduce the function 

C%(r) = m lk ffi ■ a k 5(r - r lk )r+r+\i) , 

POO 

M 2 C ? = / C%(r)dr . (11) 
Jo 

This definition is in analogy to the related function C 0u (r) first introduced in Ref. ^ anc l discussed 
in detail later. Note that while the matrix elements M 2u and get contributions only from the 1 + 
intermediate states, the function gets contributions from all intermediate multipoles. This is the 
consequence of the 5 function in the definition of C^(r). When expanded, all multipoles contribute. 
Naturally, when integrated over r only the contributions from the 1 + are nonvanishing. Examples of the 
function C 2 "(r) are shown in Fig. [3j (Similar figure appears in Ref. E7}. ) 



One can see that in all the cases shown in Fig. [3] the function C^(r) consists of a positive peak at 
r ~ 1 fm and a negative tail starting at r ~ 2 fm. It turns out that the areas under the positive peak and 
the negative tail are nearly equal, resulting in considerable uncertainty in M 2 ". 

Theoretical evaluation of the matrix elements M 2u and respectively application of the corre- 

sponding rate equation ([2]), involves another problem, namely whether one should, or should not, apply 
the experience with the ordinary (5 decay and use the concept of quenching of the GT strength <23j yj e 
shall discuss that issue in more detail later. 

4. Operator of the Of/3/3 decay 

The Ov decay rate associated with the nonvanishing value of m v is of the general form 

u 0v = 2nY> spin \R Qv \ 2 5{E el + E e2 + E f - M i )d 3 p el d 3 Pe2 , (12) 

where Ej is the energy of the final nucleus and R$ v is the transition amplitude including both the lepton 
and nuclear parts. 
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After substitution for the neutrino propagator and integration over the virtual neutrino momentum, 
the lepton amplitude acquires the form 

f d^q e~ iq ( x ~ yS) 1 1 r 

~ iS ik J ( 27r )4 q 2 _ m 2 e ( a )7pg(l - J&Wln + mj)-(l - "f 5 h*e (y) . (13) 

From the commutation properties of the gamma matrices it then follows that the decay amplitude for 
purely left-handed lepton currents is proportional to the neutrino Majorana mass rrij. Integration over the 
virtual neutrino energy leads to the replacement of the propagator (q 2 — m|) _1 by the residue ir/ooj with 
Uj = (g 2 + m J 2 ) 1 / 2 . 

Finally, the integration over the space part dq leads to an expression for the "neutrino potential" that 
appears in the corresponding nuclear transition operator, 

B( r, E m) = * / = » r *-?^L = » r . (!4) 

2TT Z J 00 00 + A m TIT Jo L0(00 + A m ) 7T Jo 9 + A m 

where the nuclear radius R = 1.2^4 1//3 fm was added as an auxiliary factor so that H becomes dimensionless. 
A corresponding 1/R 2 compensates for this auxiliary quantity in the phase space formula. The weak 
dependence on the excitation energy of the virtual intermediate odd-odd nucleus appears in A m = E m — 
Ei + E e = E m - {Mi - M f )/2. 

The momentum of the virtual neutrino is determined by the uncertainty relation q ~ 1/r, where r < R 
is a typical spacing between two nucleons. We will show later that in fact the relevant values of r are only 
r < 2-3 fm, so that the momentum transfer q ~ 100-200 MeV. For the light neutrinos the neutrino mass 
rrij can then be safely neglected in the potential H{r). (Obviously, for heavy neutrinos, with masses Mj S> 
1 GeV a different procedure is necessary.) Also, given the large value of q the dependence on the difference 
of nuclear energies E m — Ei is expected to be rather weak and the summation of the intermediate states 
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can be performed in closure for convenience. This approximation H(r,E m ) ~ H(r,E) is, in fact, typically 
used in the evaluation of the M 0u , where E ~ 5 — lOMeV is the characteristic nuclear excitation energy. 

It is worthwhile to test the validity of this approximation. Such test can be conveniently performed 
within the QRPA, where the sum over the intermediate states can be easily explicitly carried out. In this 
context one can ask two questions: How good is the closure approximation? And what is the value of the 
corresponding average energy? In Fig. [J] we illustrate the answers to these questions (see also ^3). The 
QRPA matrix elements evaluated by explicitly summing over the virtual intermediate states quoted in the 
caption can be compared with the curves obtained by replacing all intermediate energies with a constant 
E, which is varied there between and 12 MeV. One can see, first of all, that the M 0l/ changes modestly, by 
less than 10% when E is varied as expected given the relative sizes of q and E and, at the same time, that 
the exact results are quite close, but somewhat larger, than the closure ones. Thus, employing the closure 
approximation is appropriate for the evaluation of M 0li even though it apparently slightly underestimates 
the M 0v values. 




4 6 

E [MeV] 



Figure 4: Matrix elements M for several nuclei evaluated within QRPA in the closure approximation as a function of the 
assumed average excitation energy. The values of M 0v obtained without the closure approximation are 5.24 ( 76 Ge, full black 
line), 2.62 ( 96 Zr, dashed red line ), 4.99 ( 100 Mo, dot-dashed green line ), and 4.07 ( 130 Te, double dot-dashed blue line). 



The neutrino potential in the Eq. (|14p was defined assuming that the nucleons are point particles. 
That is not true, however, and thus it is necessary to include a corresponding correction in the definition 
of H (r, E). It is customary to approximate this correction in the form of the dipole type form factor 

Jfns = — — , (15) 



1 + 



M 2 , 



with Ma = 1-09 GeV. Varying Ma in the interval 1.0-1.2 GeV makes little difference. 

In addition, while the weak current of quarks has the simple V — A structure, the weak nucleon current 
contains additional terms, since nucleons are complicated composite objects, 



gv{q 2 )i^ - i9M(q 2 



a 



2m. 



-qu - 3a(<? 2 )7 m 75 + gp(q 2 )q^i5 



(16) 



representing the vector, weak magnetism, axial vector and induced pseudoscalar. In the zero momentum 
transfer limit the values of the corresponding form factors are well known. The induced pseudoscalar form 
factor is usually evaluated using the partially conserved axial- vector current hypothesis 

gA{q 2 ) x 2m p 



gp(q 



q 2 + m\ 



(17) 
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M = 1 _ I «U + \ *) + 1 Okz^W . (18) 



and the weak magnetism form factor is simply proportional to gv(q 2 ), 9m{q 2 ) = — ^n)gv(q 2 )- 

Taking these "higher order terms" into account and using the nonrelativistic limit for the nucleons, 
one arrives at the corresponding correction term gnoTil 2 ) m the main, GT part of the neutrino potential 
(for details see Ref. 1221) 

i 2 , 2/WA 2 

3 cp + m% ' 3 yep +m% ) 3\ gA(q) J ^mf, 

Finally, the neutrino potential HGT{ r ,E) governing the Gamow- Teller part of the matrix element M 0u 

with these correction factors included is of the form 

of? f°° a 

H GT (r,E 0u ) = — jo ( q r) — Vfl, NS tf)g%l T tf)d q , (19) 
TT JO q + Eq v 

When the finite nucleon size, higher order terms are neglected, and Eq u = is assumed, the potential 
has Coulomb-like shape R/r. The full potential, Eq. (|19p. however, is finite at r = 0, H{r — > 0, Eq v = 0) = 
5M^i?/16. Including the higher order currents and finite E in Eq. f|19f) increases the value of H(r = 0) by 
~30%. The shape of the neutrino potential is shown in Fig. [5j 
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Figure 5: The potential Hgt{t,E). Different approximate forms, as well as the exact one, are shown. 
In the full Qvf3f3 operator there are several potentials, each with its own spin structure 



O 0u = S iife r+r+ \-H F (r, £ ,)|r + H Gt(t, E 0v )<Ti ■ a k + H T (r, E Qu )S tk j , (20) 

where in the momentum representation = 3<7j ■ qa k ■ q — <7j • a k and g^oT^I 2 ) = 1 ; 9hot(Q 2 ) = 

o xi^'i — n ( ^1 m -2 ) ■ Also, in /fns for the Fermi part one needs to replace Ma by My ~ 0.85 GeV, and 

in the tensor part <7#ot(<7 2 ) the spherical Bessel function jo(qr) must be replaced by j2{qr). 

To finish this section, lets briefly mention the analysis of the 0v(3(3 decay assuming the existence of 
the right-handed weak currents. Such phenomenological approach is often quoted in the literature, even 
though it is not at all clear that there exist a corresponding realistic particle physics model giving the 
Oz^/3/3-decay rates competitive with the standard light left-handed Majorana neutrino exchange, like in Eq. 
Ol). The assumed hamiltonian is 



H 



G F cos 9 C 2n 
11 = ^-^= 1 



j Lin J Li + 3Ri^Ri 



(21) 
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where the lepton currents are 

j£ i = e(l- l5 )N iL ; f Rl = e(l + l5 )N m , (22) 
with the hypothetical heavy neutrinos N^m) with n mass eigenstates. The nuclear currents are 

J£l = c/eiJf; fg = v*[Xjg+T,J!*}, (23) 

where U e i is the standard light neutrino mixing matrix and V e i is its analog for the heavy neutrinos. 
Parameters A and ry characterize the strength of the right-right and right-left interactions. The formulae 
for th e corresponding Oz^/3/3-decay rate have been evaluated in many papers, here I quote just two of them 
j n those references one can find the complete expressions for all necessary neutrino potentials and 
the corresponding spin and q dependencies. 



5. Decays mediated by the heavy particle exchange 

As stressed above, when heavy particles of any kind mediate the Ov/3/3 decay, we are dealing with a 
six fermion vertex, representing extremely short range operator. Vergados ^3 was presumably the first 
author to describe how the issue of suppression, due to the short range nucleon-nucleon repulsion, can be 
overcome. 



w~ 



(a) 



(b) 



Figure 6: Schematic graph indicating the 0i//3/3 decay mediated by the heavy neutrino exchange vm in panel (a). In panel (b) 
the possible Ou/3/3 decay through the exchange of other heavy particles, in this case two squarks and a gluino in RPV SUSY. 



If one could treat nucleons as pointlike particles and assume that the heavy neutrino mass will be 
substantially larger than the corresponding momentum transfer q (say, rriheavy ^> 1 GeV) the neutrino 
potential would contain exp(-mr) where m = min(mheavy, Myy). The corresponding nuclear matrix element 
would be heavily suppressed due to the extremely short range of this potential. 

However, nucleons are not pointlike particles. Instead, their finite size could be parametrized through 
the dipole type form factor Jfns m Eq. (|15|) . In that case the neutrino potential will be of the form 



n heavyu [r,m A ) ^2 J (27r)3 I M 2 +( p I 48 



l + M A r + ^(M a r) 



(24) 



Such potentia l will have the range 1/Ma and will be much less affected by the short range correlations 
(see e.g. Ref. ESEJ^ 'pj^g disadvantage of the form factor modeling is that the error introduced by such 
approximation is very difficult to estimate. 
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In general the heavy particle exchange is characterized by some energy scale App (see Eq. ([5]), that is 
much larger than any hadronic scale Ah ~ 1 GeV that would enter the problem. In that case one could, 
instead of the form factor approach indicated above, use the prescriptions of effective field theory and 
classify the contributions in powers of small quantities q/Ajj, qj App and Ah/ App as in Ref. ^2}. 

The hadronic vertices appearing in the corresponding Lagrangian will be of the type NNNNee, NNiree 
and mree as illustrated in Fig. [7J They stem from quark-lepton operators having different transformation 
properties under parity and chiral SU(2). As such, they will contribute to different orders in the q/An 
expansion. The vertices involving pions are longer range. They have been analyzed in the form factor 
approach in Ref. but the EFT allows mor e sy stematic approach because of the separation of scales 
q < Ah <C App. It was noted already in Ref. ^^it that the nuclear matrix elements associated with the 
long range pionic effects within the RPV SUSY scenarios can be dominant. But that is, in fact, a more 
general result. The pionic effects can be substantially larger than those obtained using the conventional 
form factor model for the short-range NNNNee process. 



Figure 7: Diagrams that contribute to Ovfifi decay at tree level. Panel (a) represents nuee contribution, panel (b) NNiree 
contribution and (c) NNNNee contribution. 
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The lepton number violating vertices in Figfjj represent nonstandard model operators, and can be char- 
acterized by some parameters K n7T , K^n^ and K^nnn- Counting powers of qj Ah, for Ah = 4irf w , f w ~ 
92 MeV, we come to the conclusion that 



Fig. [Th, ~ K^p 2 , Fig. [7b, ~ K NN7r p \ Fig. [7b ~ K NNNN p 



o 



(25) 



Thus, the long range Oz^/3/3-decay operators iriree and to a lesser degree NNiree are enhanced, relative to 
the short range operator NNNNee in FigJTJc). 

Clearly, the most important operator is the one corresponding to the iriree vertex. It was shown in 
Ref. ^ that the corresponding Lagrangian can be expressed as 



L (0) 



A 



{vr-vr-e^! + /3 27 5 )e c + vr+vr+e^A - /3 2 7 5 )e} 



(26) 



where /?i 2 are dimensionless parameters that need be evaluated in any concrete particle physics model. 

Transforming the above Lagrangian into the hadron-lepton system, taking the nonrelativistic limit and 
Fourier transforming to coordinate space yield an expression for the transition matrix element 



Mo 

Sin 



R 



(27) 



Here t%« is the distance between the zth and jth neutrons, r = r/r and X{j — TijTfl^. The form-factors Fi 



(28) 



and F2 were first introduced in Ref. ^ 

Fi(x) = (x-2)e~ x , F 2 (x) = (x+l)e- 
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The Ovf3f3 halflife is, as expected, inversely proportional to the square of the scale and depends on 
the empirical parameters f3\ 2 , 



1 fir 2 n 4 A 4 H A r E 00- m e 



^[(ft 2 + f3 2 2 ) P iE lP2 E 2 - {Pi - (3 2 ) Pl p 2 m 2 e ]\M \ 2 , (29) 

As mentioned earlier, concrete application of the pion exchange mechanism in Ref. ^3 suggests the 
dominance of the iriree over the short range nucleon only vertex by a factor of 10 - 30 at the level of the 
nuclear matrix element. This is, to some extent, so far largely unexplored issue. The application of the 
form factor method, obviously insufficient, is more or less a norm even today (see e.g. the recent paper 



6. Short range correlations 




distance r (fm) 



Figure 8: Plots of the squares of the Jastrow-like correlation functions described in the text. 

When evaluating the Qv(3(3 nuclear matrix element in the closure approximation it is necessary, ulti- 
mately, to consider the matrix element of a two-body transition operator which connects a state with two 
neutrons in a specific mean field single particle states with two protons again in some specific final states. 
These states are angular momentum and possibly isospin coupled single particle states. In essentially all 
cases they do not explicitly include the effect of the short range nucleon- nucleon repulsion. 

Until recently, it was customary to simulate the effect of such short range correlations by including in 
the operator (or equivalently in the two body wave-function) a phenomenological Jastrow-like function 

fjastro W (r) = 1 - ce- ar2 (1 - br 2 ) . (30) 

In other words, 

M 0u = Yj k<m M km M km -)■ f(r km )M km f(r km ) (31) 

In early works the parameters a, b, c, determined for a different purpose in Ref. were used (c = 1.0, a = 
1.1 fm -2 and b =0.68 fm -2 ). This is the so called Miller-Spencer parametrization. Since in the evaluation 
of the nuclear matrix element such function appears twice, we plot in Fig. [8] its square. Application of 
this prescription resulted typically in the reduction of the corresponding 0(v/3/3 nuclear matrix element by 
~30% when compared to the result when fj as trow(r) = 1 is used. 

Recent works ha ve q uestioned the adequacy of this prescription. One of the approaches proposed 
relatively recently in i s based on the Unitary Correlation Operator Method (UCOM), see Ref. S3. 
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This approach describes not only short-range, but also central and tensor correlations explicitly by means 
of a unitary transformation. Applied to a realistic NN interaction, the method produces a correlated 
interaction, which can be used as a universal effective interaction. In the case of the O^/3/3-decay calculation 
the correlated two-nucleon wave function was taken as 



l*7> = C r \9j). (32) 

Here, C r is the unitary correlation operator describing the short-range correlations. The explicit form 
of C r is given in The UCOM-corrected nuclea r matr ix elements are significantly less reduced when 
compared with those calculated with Jastrow SRC j n p rac tise, the UCOM prescription requires 

that the operator r a b be replaced by a shifted R + (r a b). As pointed out in Ref. ^D, care must be taken to 
apply the p rescrip tion consistently, not just simply by replacing the neutrino potential H(r) — > H(R + (r)). 

In Refs. t he effect of short range correlations was computed within well defined Brueckner-based 

approximation scheme. In particular, in the coupled cluster method was used, based on the realistic 
NN interactions, CD-Bonn and Argonne V18. In a good approximation the effect of correlations can 
be again approximated by the Jastrow- like functions (see Eq. (|30p ) with parameters fitted in Ref. ^2}, 
a = 1.59 fm -2 , b = 1.45 fm -2 and c=0.92 for the case of the Argonne potential and a = 1.52 fm -2 , b = 1.88 
fm~ 2 and c = 0.46 for the case of the Argonne potential. Both of these functions are displayed in Fig. 
It is obvious that the application of these functions leads to the much less reduction of the M 0u , in 
agreement with the UCOM method and with the conclusion of Ref. S3 than the Miller-Spencer approach 
used previously. In fact, provided the nucleoli finite size form factor are properly taken into account, the 
application of the UCOM procedure or the prescripition of Ref. ^3) leads to the M 0u values essentially 
unchanged in comparison to the omission of the short range correlation correction altogether. Thus, the 
M Gu are about 30% larger than after the application of the traditional Miller-Spencer parametrization. 
There is a general consensus that the matter of treatment of the short range correlation is satisfactorily 
resolved at the present time. 



7. nuclear matrix elements: Basics 



In this section the basic issues involved in the numerical evaluation of the Oz^/3/3 nuclear matrix elements 
will be described. We concentrate here on the standard light Majorana neutrino exchange mechanism, so 
that the relation between the decay rate and the nuclear matrix elements is governed by Eq. ([2]). First, 
there are few reasonable assumptions common to all (or essentially all) methods: 

• In the + — > + transitions the outgoing leptons are in the S\/2 (or more precisely Dirac's k = — 1) 
state. That assumption is used in the calculation of the phase space factor. 

• The hadronic currents are treated in the nonrelativistic impulse approximation. 

• In most, but not all, applications the closure approximation is used, i.e. the matrix element of a 
two-body operator is considered (see Fig. 0|). 

Since the decay involves the transformation of two bound neutrons into two bound protons, it is 
necessary, first of all, to choose the proper mean field in which the nucleons are bound. That field could 
be spherical, but it can be also deformed. In the case of deformed mean field the corresponding intrinsic 
states have no definite angular momentum; projection into states with definite value of J n is sometimes 
performed in that case. 

Once the mean field has been specified, the set of single particle states \ipi) is obtained. The next step 
is to decide which of these single particle states are fully occupied or totally empty in both the initial and 
final nuclei. Such states, obviously, cannot participate in the decay; the fully occupied states form the 
inert core. The remaining states, i.e. the single particle states that are partially occupied or where the 
occupancies in the initial and final nuclei are different, form the valence space. To obtain the resulting 
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many-particle states one has to take into account the residual interaction among the nucleons constrained 
in the valence space. 

While such division into occupied, valence, and empty states appears to be reasonable and well defined, 
in practice real nuclei have diffuse Fermi levels and it is not a priori clear how to properly decide to which 
category a given single particle state belongs. Moreover, the effect of the core and empty states should 
be included, in principle, in the renormalization of the nuclear hamiltonian and in the definition of the 
effective operators. See Ref. anc i the contribution of J. Engel to this focus issue for explanation and 
references; there are not many explicit applications of this procedure to the problem at hand, so far. 
Ideally, moreover, the mean field should be determined selfconsistently, so that the same hamiltonian, or 
nucleon-nucleon interaction, is used in its definition and in the treatment of the interaction of valence 
nucleons. Note that for ease of computation the single particle wave functions are typically taken to be 
the eigenfunctions of the harmonic oscillator potential, or superpositions of such functions. 

Variety of methods has been used for evaluation of the Ov/3/3 nuclear matrix elements. They differ 
in their choice of the valence space, interaction hamiltonian and the ways the corresponding equations of 
motion are solved. Some of these methods are described in detail in separate contributions to this focus 
issue. Here I wish to characterize each method very briefly and show its typical output. Let me stress that 
an exact "ab initio" , i.e. without approximations, calculation of M 0u for the candidate nuclei is impossible 
at the present time. 

The nuclear shell model (NSM) is, in principle, the method that seems to be well suited for this task. In 
it, the valence space consists of just few single particle states near the Fermi level (usually one main shell). 
With interaction that is based on the realistic nucleon-nucleon force, but renormalized slightly to describe 
better masses, energies and transitions in real nuclei, all possible configurations of the valence nucleons are 
included in the calculation. The resulting states have not only the correct number of protons and neutrons, 
but also all relevant quantum numbers (angular momentum, isospin, etc). For most nuclei of interest ( 48 Ca 
is an exception) the valence space, however, does not include enough states to fulfill the Ikeda sum rule 
(see Eq. (|10p ). hence full description of the /3 strength functions is not possible. However, NSM is well 
tested, since it is capable to describe quite well the spectroscopy of low lying states in both initial and final 
nuclei. In the following figures 191101 and [11] the NSM results are denoted by the blue squares. 

The 2u/3/3 decay matrix elements M 2u for several nuclei in Table [T] are reasonably well described in 
the NSM, see Ref. SI (ioo Mo being a notable exception). However, to achieve this task, it was necessary 
to apply quenching factors that, for nuclei heavier than 48 Ca, are considerably smaller than in the lighter 
nuclei where the valence space contains the full oscillator shell. Note that no quenching is applied to the 
results shown in Figures 191101 and [TT1 I will describe the issue of quenching of the weak nucleon current 
operators in Section,, 

The quasiparticle random phase approximation (QRPA) and its renormalized version (RQRPA) is 
another method often used in the evaluation of M 0u . In it, the valence space is not restricted and contains at 
least two full oscillator shells, often more than that. On the other hand, only selected simple configurations 
of the valence nucleons are used. The basis states have broken symmetries in which particle numbers, 
isospin, and possibly angular momentum are not good quantum numbers but conserved only on average. 
After the equations of motion are solved, some of the symmetries are partially restored. The RQRPA 
partially restores the Pauli principle violation in the resulting states. 

The procedure consist of several steps. In the first one the like particle pairing interaction is taken 
into account, using the BCS procedure. Then, the neutron-proton interaction is used in the equations 
of motion, resulting in states that contain two quasiparticle and two quasi hole configurations and their 
iterations. Usually, the realistic G-matrix based interaction is used, but block renormalized with common 
renormalization factors. In particular, all particle-particle channel interaction matrix elements are typically 
scaled by the factor g pp with nominal, unrenormalized value g pp = 1. In many applications, including the 
one used in Figsl9| [TOl and [TT| the parameter g pp is determined such that the experimental matrix element 
M 2u is correctly reproduced. This procedure has been first suggested in Ref , resulting in <20% change 
in g pp . In this case, the method obviously is not predicting the magnitude of M 2u , instead it uses its 
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Figure 9: Dimensionless 0^/3/3 nuclear matrix elements for selected nuclei evaluated using a variety of indicated methods. For 
references see text. 

value in the fit. However, as shown inS3 ; the procedure removes the dependence of M 0u on the size of 
the valence space, making its predicted value essentially constant. While most early QRPA and RQRPA 
calculations assume spherical nuclear shape, recent Ref. S3 extension allows deformed nuclear shape. In 
Figs. [U [10] and HU the RQRPA results are denoted by red circles. 

The IBM-2 method uses the microscopic interacting boson model to evaluate M 0u S3. i n IBM-2 
one begins with correlated S and D pairs of identical nucleons and includes the effect of deformation 
through the bosonic neutron-proton quadrupole interaction. The method describes well the low lying 
states, the electromagnetic transitions between them and the two-nucleon transition rates in spherical and 
strongly deformed nuclei. Even though the method was originally considered as an approximation of the 
nuclear shell model, the resulting M 0l/ (purple upward pointing triangles in Figs. l9l 1101 and lllj) are rather 
surprisingly close to the RQRPA results and noticeably larger than the NSM ones. At the present time it 
is not possible to evaluate the M 2u matrix elements within the IBM-2 method. 

The projected Hartree-Fock-Bogoljubov framework SHI (PHFB) uses angular momentum projected 
wave functions based on the hamiltonian with pairing and quadrupole-quadrupole interaction. Again, the 
method is not capable in describing the intermediate odd-odd nuclei and thus neither the M 2u matrix 
elements. The resulting M 0u matrix elements for the heavier nuclei (including the deformed 150 Nd) are 
denoted with brown downward pointing triangles in Figs. [10] and [UJ This method, like IBM-2, does not 
explicitly contain the isoscalar neutron-proton interaction, crucially important in QRPA. Yet its resulting 
Al 0u are in a reasonable agreement with those from RQRPA. 

In Ref. S3 the generating coordinate method (GCM) was used in conjunction with the particle number 
and angular momentum projection for both the initial and final nuclei. Large single particle space was 
used (11 shells) with the well established Gogny D1S energy density functional. The initial and final 
many-body wave functions are represented as combinations of the particle number N, Z projected I = + 
axially symmetric states with different intrinsic deformations. The lowest state is found by solving the Hill- 
Wheeler- Griffin eigenvalue equation. The resulting M 0u matrix elements are denoted with green diamonds 
in Figs. 19)101 and 111! Again, like in PHFB and IBM-2 methods, the treatment of the odd-odd nuclei, and 
thus also of the M 2u matrix elements is impossible in GCM. 

Looking at the Fig. [9] one can see that, common to all five displayed methods the predicted M 0u nuclear 
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Figure 10: Halflives T^n in units of 10 27 years corresponding to the effective neutrino mass < mpp > = 20 meV and the 



matrix elements shown in Fig. [9] The lifetime T1/2 scales as < mpp >" 
is not shown. 
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Figure 11: The effective neutrino mass < mpp > in units of meV evaluated assuming that the halfiife of all nuclei is 10 27 years 
and for the matrix elements shown in Fig. [9] Note that < 771,3,3 > scales as T~^J 2 . 
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matrix elements vary relatively smoothly, with the mass number A, unlike the experimentally determined 
M 2u matrix elements displayed in Fig. [2] that vary strongly with A. The RQRPA, IBM-2, PHFB and 
GCM method are in a crude agreement with each other, and predict slow decrease of M 0u with increasing 
A. The M 0l/ evaluated in NSM are essentially constant with A and noticeably smaller than those from 
the other methods, particularly in the lighter nuclei. In Figs. [TO] and [TT] I show the halflives Tu 2 for the 
fixed effective neutrino Majorana mass {mpp) and {mpp) for fixed TV^. These figures should be helpful for 
comparing the experimental results and/or sensitivities for different nuclei. 

The question that is, for obvious reasons, often asked is what is the error or uncertainty of individual 
nuclear matrix elements? It is notoriosly difficult to estimate the error of a theoretical result. One 
possibility is just to use, without selection, all calculated values and take as the error their spread. That 
was used in the provocative Ref. ED. Clearly that is not the correct approach, even though the quoted 
paper served as appropriate warning to the nuclear structure community. 

Another possibility is to use as a measure of uncertainty the spread of selected careful calculations, like 
those in Fig. [9] That also has problems. Different methods use different approximations, and it is unlikely 
that all of them are equally important. And it is difficult to decide which of them is more realistic and 
which is less realistic. 

Within a given nuclear structure method it is possible to assign a measure of error by considering the 
uncertainties in the input parameters. That was done for QRPA e.g. in Refs. | 45 | 40 t w here the size of 
the single particle basis, whether QRPA or RQRPA was used, and whether g^ was quenched or not is 
used for such estimate. Clearly that error does not include the possible systematic uncetainty of the basic 
method itself, but serves as a useful estimate of the irreducible uncertainty. When considering the ratios 
of the matrix elements for different nuclei, the question of correlation of the corresponding errors arises. 
As shown in Ref. ^ at least in the case of QRPA, the errors are highly correlated. 

8. Qi//3/3 nuclear matrix elements: Physics considerations 

Double beta decay in both 2u and Of modes can exist because even-even nuclei are more bound than 
the neighboring odd-odd nuclei. This extra binding is a consequence of pairing between like nucleons. In 
nonmagic systems neutrons and/or protons form + pairs and the corresponding Fermi level becomes diffuse 
over the region with the characteristic size ~ pairing gap A. This opens more possibilities for nn — > pp 
transitions. The calculated matrix elements M 0u increase when the gap A increases. The (unrealistic in 
real nuclei) situation of pure paired nuclear system (only seniority states) would have very large M 0u . 

However, in real nuclei opposite tendencies exist. Real nuclei have admixtures of the "broken pair" 
states, or in the shell model language, states with higher seniority. These states are present because other 
parts of the nucleon interaction exist, in particular the neutron-proton force. It is illustrative to characterize 
such states by the angular momentum J of the neutron pair that is in the j3j3 decay transformed into the 
proton pair with the same J . In Fig. [12] the corresponding competition is illustrated, in both NSM and 
QRPA. While the positive pairing parts are very similar (since the same single particle spaces are used), 
the negative J ^ parts are only qualitatively similar. Nevertheless the severe cancellation between these 
two tendencies, and between the corresponding components of the residual interaction, is present in both 
methods. 

The situation shown in Fig. [12] is somewhat unrealistic as far as the QRPA is concerned; the corre- 
sponding single particle space is just too small. When it is enlarged to contain two full oscillator shells, the 
competition becomes even more severe. Both the J = and J ^ grow, while their difference remains 
more or less unchanged. This situation, illustrated in Fig[13j appears to be is universal. One can see 
that, due to this cancellation, accurate evaluation of M is difficult; small uncertainty in either of the two 
(positive and negative) components is strongly enhanced in the difference. 

The Ozy/3/3 decay operators depend on the internucleon distance ryi due to the neutrino potential 
H(r, E). Obviously, the range of ryi is restricted from above by r\2 < 2R. From the form of H (r) ~ R/r 
one could, naively, expect that the characteristic value of r\i is the typical distance between nucleons in 
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Figure 12: Contributions of different angular moments J of the two decaying neutrons to Mgy in 82 Se (upper panel) and 
130 Te (lower panel). The results of NSM (red) and QRPA (blue) are compared. Both calculations use identical limited single 
particle spaces. The values of Mg T are 1.32(QRPA) and 2.06(NSM) for 82 Se and 1.05(QRPA) and 1.98(NSM) for 130 Te. 




Figure 13: Contributions of different angular moments J of the two decaying neutrons to Mq" t in 76 Ge (blue), 100 Mo (red) 
and 130 Te (green) within QRPA with proper single particle space (two full oscillator shells). The M°" values are 3.98 ( 76 Ge), 
2.74 ( 100 Mo and 2.67 ( 130 Te. 
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the nucleus, namely f\i ~ R. However, that is not true, in reality only much smaller values of r\2 < 2-3 fm 
or equivalently larger values of the momentum transfer q are relevant. That was first demonstrated within 
QRPA in Ref . SDJ ; very similar result was obtained also in NSM EZJ _ 
To see how that conclusion is obtained, define a function C(r) 

CcAr) = • d k T+r+S{r - r lk )H(r lk ,E)\i) , (33) 

Obviously, this function is normalized by 

M°Jt = / C°J T (r)dr , (34) 
J o 

and is related to the function C^(r) (see Eq. ([11]) ) by 

C°J T (r) = H(r,E)xC 2 J(r) , (35) 

which is valid for any shape of the neutrino potential. 

Examples of the functions C 0u (r) are shown in Fig. [TUfor three representative nuclei. As the lower 
panel demonstrates, the cancellation between the "pairing" (J = 0) and "broken pairs" {J ^ 0) is 
essentially complete for ryi > 2-3 fm (similar figure appears in®-). That cancellation is so complete 
when the particle-particle renormalization constant g pp is chosen in the usual way, i.e. so that the 2v(3(3 
lifetime is correctly reproduced. For other values of g pp the cancellation between J = and J ^ is 
less perfect. Analogous conclusions can be obtained in an exactly solvable model based on the algebra 
SO{5) x 50(5). 

The relation between C 0u (r) and C^(r) in Eq. (|35p is intriguing. Can it be used for the evaluation 
of M 0u ? There are two problems with that possibility. While M 2u can be deduced from the measured 
lifetime of the 2vj3j3 decay, the closure matrix element cannot be deduced that way. Moreover, even if it 
would be possible to somehow determine accurately, its magnitude is insufficient to determine the 
whole function C^(r). 



9. Quenching of the weak axial current 

It is well known that experimental Gamow- Teller /3-decay transitions to individual final states are no- 
ticeably weaker than the theory predicts. That phenomenon is known as the axial current matrix elements 
quenching. The /3-strength functions can be studied also with the charge exchange nuclear reactions and 
a similar effect is observed as well. Thus, in order to descr ibe the m atrix elements of the operator <tt, 
the empirical rule (o't) 2 ^ ~ 0.6(<7T)^ iorfei is usually used (see^2Ell55j-j_ Si nce these operators accompany 
weak axial current, it is convenient to account for such quenching by using an effective coupling constant 
g e j[ ~ 1-0 instead of the true value gA = 1.269. 

The evidence for quenching is restricted so far to the Gamow- Teller operator <tt and relatively low-lying 
final states. It is not known whether the other multipole operators associated with the weak axial current 
should be quenched as well. In fact, the analysis of the muon capture rates in ESJ an( j Q £ ^ e unique second 
forbidden /3 decays in ^3 suggests that quenching is not needed in those cases. 

The quenching of the axial current matrix elements, or simply of gA, is believed to be caused by nucleon 



correlations 1 ^'. The alternative explanation, screening of the GT operato r by the A-hole pairs ^3 is, 
these days, not considered as the likely explanation. Very recently, in Ref. the quenching has been 
associated with the two-body currents. 

Since the 2vj3j3 decay involves only the GT operators and relatively low-lying intermediate states, one 
could expect that the quenching is involved in that case. Indeed, as already mentioned, in the nuclear 
shell model the agreement with experimental decay rate is achieved only with quenching (ar) 2 ^ ~ (0.2 — 
0.36)(aT) 2 model , with similar quenching required to reproduce the measured /3 decay rates and the f3 strength 
functions. 
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Figure 14: Dependence on ryz for 76 Ge (full line), 100 Mo (dot-dashed line) and 130 Te (dashed line) evaluated in QRPA. The 
upper panel shows the full matrix element, and the lower panel shows separately the "pairing" [J = 0) and "broken pair" 
( J 7^ 0) contributions. The integrated M°" are 5.4, 4.5, and 4.1 for these three nuclei. 



In QRPA and RQRPA the usual way of the g pp renormalization using the experimental M means 
that the effect of quenching cannot be deduced from the comparison of the calculated and empirical M 2u . 
In the other discussed methods (IBM-2, PHFB, GCM) it is not possible at the present time to calculate 
M 2v ', so the question of quenching for the 2^/3/3 is meaningless in that case. 

However, whether quenching should be applied to the 0vf3f3 matrix elements M® u is an important issue 
and a source of noticeable uncertainty. Within Q RPA and RQRPA the effect of quenching on the Oz^/3/3 
decay mode has been considered in Refs. I 27 | 43 | 61) _ The calculations there uses a modified definition 

M'°» = (Jj^j M°"(gf), (36) 

and thus the standard phase space factor can be used independently of whether quenching is included or 
not. It is concluded there that the matrix elements M' 0u are reduced by ~20-30% when g e P = 1.0 is 
assumed compared with gA = 1.269. Thus, in that case the predicted 0v decay rates are affected by the 
possible quenching less than naive expectation based om the ratio [<?^ /<m] 4 might suggest. 

Somewhat similar, but from the point of view of physics involved quite different reasons was reached 
in Ref. when the effect of two-body current is included in the effective field theory. 

The problem of quenching of the M 0u matrix elements remains one of the main sources of uncertainty 
in their true value. The relative change with A and Z, as depicted in Fig. [9] will be, however, affected only 
little by the inclusion of the quenching phenomenon. 

10. Conclusions 

In the last few years the problem of evaluation of the nuclear matrix elements M 0l/ of 0u/3/3 decay has 
entered a new era. As described here, there are now five or more seemingly quite different methods of 
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solving it; it became one of the very active subfields of nuclear structure theory. The results are not yet 
quite convergent, but at least they suggest that the calculated values are fairly insensitive (within a factor 
of about two) to the broad range of approximations made. Moreover, all calculations agree that the values 
of M Qu do not change abruptly from one candidate nucleus to another one. Thus, if the 0uf3f3 decay were 
observed in one nucleus, one can with some confidence predict its lifetime in the other candidate nuclei, 
increasing the chances that a reliable and confirmed result is obtained. Of course, the fact that the nuclear 
shell model consistently predicts M 0l/ values that are noticeably smaller than in the other methods remains. 
Is it the inclusion of complicated configurations (states with high seniority) that causes the reduction of 
M 0l> in NSM, or is it the inclusion of larger single-particle spaces in most other methods? The answer to 
this question is obviously needed. 

Moreover, there are several more general issues that deserve closer attention of nuclear theorists. One of 
them is the problem of quenching of the axial weak current matrix elements. That effect is well established 
in ordinary allowed beta decays, and in the evaluation of the beta strength functions involving low-lying 
nuclear states. Should an analogous effect be applied to the theoretical evaluation of M 0u l This is, as yet, 
an open question. Its solution could affect the true values of M 0l/ appreciably, by about 30% according to 
our estimates. 

Problem of quenching is just a part of the determination of effective operators. How does one con- 
sistently include the renormalization caused by using only a subset of the general Hilbert space in the 
calculational framework? Renormalization of the gA, i-e. inclusion of quenching, is an example of such 
effect. The problem of taking into account the high momentum part of the nucleon interaction, causing 
the short range repulsion, belongs to that category as well. There is a consensus now, that the recent 
developments (see Section ^ point us in the right direction. More work in that direction would be clearly 
beneficial. 

Finally, most attention up till now was concentrated on the so-called standard scenario, according to 
which the Ov/3/3 decay would be caused by the exchange of the light Majorana neutrinos that interact 
through the canonical left-handed weak currents. Much less attention has been given to the possible 
mechanism involving heavy, ~ TeV, particle exchange and thus extremely short range effects. Is it really 
generally true that in these cases the pionic effects, described in Section „ dominate? If that is the case 
much more detailed evaluations of the M 0l/ in such cases is clearly needed. In particular, in the next few 
years much progress, one hopes, in the exploration of the TeV mass range particles will be achieved at LHC. 
If some of the suggested particle physics models (Left-Right symmetry, R-parity violating supersymmetry) 
might find support at LHC, more work on the corresponding M 0u matrix elements will be clearly needed. 
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Appendix: Possible resonances in OuECEC decays 

The two electron capture decay without neutrino emission requires a special comment. Clearly, when 
the initial and final states have different energies, the process cannot proceed since energy is not conserved. 
The radiative process, with bremsstrahlung photon emission, however, can proceed and its rate, unlike all 
the other neutrinoless processes, increases with decreasing Q value i§2& . (However, the estimated decay 
rates are quite small and lifetimes long). In the extreme case of essentially perfect degeneracy the photon 
emission is not needed, and a resonance enhancement can occur ^3 . 

In the case of resonance the initial state is the atom (Z, ^4), stable against ordinary f3 decay. The final 
state is the ion (Z — 2, ^4) with electron vacancies H, H' and, in general, with the nucleus in some excited 
state of energy E*. The resonance occurs if the final energy E = E* + Eh + E^ 1 is close to the decay Q 
value, i.e. the difference of the initial and final atomic masses, and a perfect resonance occurs when Q — E 
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is less than the width of the final state which is dominated by the electron hole widths Th, Tjj 1 - The decay 
rate near resonance is given by the Breit-Wigner type formula 

1_ (AA/) 2 „ 



t (Q - B)' + 1-2/4 



where AM is the matrix element of weak interaction between the two degenerate atomic states. 

The states of definite energy, the eigenstates of the total hamiltonian, are superpositions of the initial 
and final states, mixed by AM. But in reality, the initial state is pure, and not a state of definite energy, 
since the final state decays essent ially immediately. 

The mixing matrix element is ^* 

AM ~ ^^{m^\m?9\M^ , (38) 

where ip(0) is the amplitude at the origin of the wave function of the captured electrons and M® u is the 
nuclear matrix element, same one as before. Clearly, if the resonance can be approached, the decay rate 
would be enhanced by the factor 4/T compared to T/(E — Q) 2 , where the width T is typically tens of eV. 
Estimates suggest that in such a case the decay lifetime for {mpp} ~ 1 eV could be of the order of 10 24-25 
years. However, chances of finding a case of a perfect (eV size) resonance when E is of order of MeV are 
not large. Indeed, in the best case found so far, in 152 GdE0, the quantity Q- E = 0.91(18) keV, still 
with the predicted halflife of only ~ 10 26 y for {mpp) = 1 eV. 
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